Trapped-ion qutrit spin molecule quantum computer 
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We present a qutrit quantum computer design using trapped ions in the presence of a magnetic 
field gradient. The magnetic field gradient induces a "spin-spin" type coupling, similar to the 
J-coupling observed in molecules, between the qutrits which allows conditional quantum logic to 
take place. We describe in some detail, how one can execute specific one and two qutrit quantum 
gates, required for universal qutrit quantum computing. 
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Quantum information and quantum computing p], 
have made huge advances both theoretically and experi- 
mentally in recent times. There are many different pro- 
posals 0, H H H |E II !, H E Et for the physical 
implementation of a quantum computer, all of which 
are specified by the physics of their qubit systems and 
the nature of the interactions between qubits. The lat- 
ter is necessary for the execution of conditional qubit 
logic, a requirement for universal qubit quantum com- 
putation. The qubit is of course the obvious unit of a 
quantum computer given our classical computer histor- 
ical dependance on binary logic. Quantum algorithms 
and protocols for quantum communication and cryptog- 
raphy have been studied extensively with the qubit as the 
information storage and transport medium. Recently d- 
level quantum systems, or qudits, have started to be con- 
sidered seriously in terms of generalizing and improving 
qubit-based quantum algorithms and protocols. Interest- 
ing results have emerged in the qutrit [d = 3) case. It 
has been shown that quantum cryptography protocols are 
more robust against eavesdropping attacks when qutrits 
are used [T^ . ll3lll4| . Quantum bit commitment and coin- 
flipping protocols are more secure with entangled qutrits 
than with qubits [l^. Indeed, it is expected that qutrit- 
based quantum information processing will be more pow- 
erful than other qudit implementations since they opti- 
mize the Hilbert space dimensionality ^t|- More spec- 
ulative is the possibility of a higher error-tolerance for 
fault-tolerant qutrit quantum computation |18|. How- 
ever, there have not been many qudit-based quantum 
computer proposals. As far as we are aware there has 
been only one qutrit quantum computer proposal using 
trapped ions which generalises the original Cirac-Zollcr 
design [l9| . 

Here we describe a modification to previous work 
pfl EH |22| , on ion trap quantum computers where now 
qutrits store the quantum information. An axial mag- 
netic field gradient is applied across an ion chain that 
allows the three hyperfine Zeeman energy levels of each 



ion, forming the qutrit, to be individually frequency ad- 
dressed. It also introduces an inter-qutrit coupling that 
facilitates conditional quantum logic between qutrits. 
Previously the operation of an ion trap quantum com- 
puter in the presence of a magnetic field gradient has 
been discussed with qubits as the unit of quantum in- 
formation. In [2fl EH], it is shown that all quantum gate 
operations, normally requiring optical irradiation, can be 
implemented using long wavelength radiation due to the 
effects of the magnetic field gradient and trapping poten- 
tial. The gradient also introduces a term in the Hamilto- 
nian that is analogous to the spin-spin coupling observed 
between nuclei in molecules in NMR. This coupling can 
be used to perform quantum logic. This idea is investi- 
gated further for ions in a linear array of microtraps 
In the only other qutrit ion trap quantum computer pro- 
posal |19| , the quantized collective vibrational motion of 
the linearly trapped ions is used as a quantum bus to 
preform quantum logic. 

In this work we consider N ions in a linear ion trap 
in the presence of a magnetic field gradient. Three un- 
equally spaced hyperfine Zeeman levels serve as our qutrit 
(Fig[TJ. We will refer throughout to the 171 Yb + ion as an 
example qutrit, particularly the F = I hyperfine Zeeman 
levels shown in Fig. [21 Our computational basis is now 
{|0),|f),|2)} written as |0) = (0,0, 1) T , |1) = (0,1, 0) T 



and |2) = (1,0, 0) T . We denote by 
frequencies resonant with the |0) <-> 
transitions for the nth ion and let A n 
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For each ion these frequencies, due 
to the spatial dependence of the magnetic field, are a 
function of the ions' positions. The internal electronic 
Hamiltonian describing the spin degrees of freedom of a 
single ion is given by H sp ^ n = \hZ n with Z n expressed 
in the above computational basis as 



Z n — 
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As in |20ll21| . the ions sit in an effectively Id harmonic os- 
cillator potential along the trap axis and feci their mutual 
Coulomb repulsion. Expanding this potential around 



2 



their equilibrium positions allows their vibrational mo- 
tion to be treated collectively and the Hamiltonian for 
the motional degrees of freedom of the ions in normal 
coordinates is 



(2) 



The local and normal coordinates are related by q — 
DQ, with D the unitary transformation that diagonalises 
the Hessian, A, of the ions potential evaluated at the 
equilibrium positions of the ions, zo yn [See refs pll Esl]]. 

The magnetic field gradient, = b, means the ions 
feel a spatially varying magnetic field, B(z) = Bq + b ■ z. 
This introduces a new term into the spin part of the 
Hamiltonian, 
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Setting M n = dZ n /dz\ Zo n for easier notation, as long as 
< ^KM n q n > is much smaller than the ground state en- 
ergy of the collective vibrational motion, then this new 
term will have a negligible effect on the normal modes 
and can be treated as a perturbation. This places a limit 
on the size of the magnetic field gradient but is no more 
stringent than other constraints discussed later. The re- 
sultant Hamiltonian is 
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We move to a rotating frame where the spin and motional 
degrees of freedom are decoupled via the unitary trans- 
formation U = cxp i J2 n M n Din^j Pq^i , yielding 
H = UHU\ with 
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where H M m = f J2i v i [y^f E„ ^nA n ) , and the po- 
sition and momentum operators have been expressed 
in terms of creation and annihilation operators in the 
usual way. Hmm can be expressed as Hmm = 



-T 

2 Z->r. 
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where 



The Hamiltonian H in (JSJ, describes iV individually ad- 
dressable qutrits coupled through Hmm, a "spin-spin" 
type interaction which we will presently show can be used 
to perform conditional quantum logic between qutrits. 
There are no extra experimental requirements compared 
to the setup proposed in pol l2l| . 

It is necessary to demonstrate how single qutrit opera- 
tions and conditional logic, the most basic ingredients of 
a universal quantum computation, can be implemented 
in this modified design. 

In order to perform single qutrit operations, we use the 
fact that the two transitions |0) <-» |1), and |1) |2), 
have different resonant frequencies allowing the opera- 
tions Uqx, and U12, 
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to be carried out. Since the |0) <-> \2) transition is for- 
bidden, any rotation between these two states requires a 
7r-pulse between |0) and |1), J7oi(f ,0), the required gate 
on |1) |2), Ui2(9,4>), followed by another 7r-pulse on 
|0) «-* |1). Also required is the unitary differential phase 
operation, 



(8) 



where this gate is a composition of "cr z " operations on the 
|0) «-» |2) and |0) <-» |1) transitions that are themselves 
compositions of U12 and Uqi operations. In particular, 
U D = (Z 02 ) p (Z 01 ) a with (Zij) p = 1 1,, I ,,.(/;. 0)//^ and 
Hij — U{j(j, f ) is the Hadamard gate. Thus the unitary 
operations given in Q allow us to generate any operation 
in SU{3) ||. 

For gates between more than one qutrit, we propose to 
use the last part of the Hamiltonian H in J^}. We now 
consider in more detail the hyperfine Zeeman levels for 
the qutrits. For intermediate magnetic field strengths 
B, such that gj^sB « A, where gj is the Lande g- 
factor, fig, is the Bohr magneton and A, is the hyperfine 
constant, the energy levels are described by the Rabi- 
Breit formula (25|. In this case dA n /dz — gj^Bb/h, and 
dS n /dz « gj/iBb/V2h, and where b is the constant gra- 
dient of the magnetic field B(z) = Bq + bz. We absorb 
the gj^sb/h factor into the definition of J nm and write 
M = dZ/d z\ Zq in the matrix form 
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(6) with now J nm = (gJ ^ b) (A 1 ) nm [13 . Expressing the 
operator M in @ in terms of the generators of SU (3) 
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we have 



where 



M = a I + a 3 A 3 + a 8 A 8 
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575' a3 = - X )> and ° 8 = 57^ 

3-\/2). Hmm is thus a two-body N-qutrit Hamiltonian 
as defined in |26j . and so, as arbitrary local unitaries 
are possible, universal quantum computation can be per- 
formed. More specifically, it has been pointed out in ^| 
that a generalized XOR gate between qutrits, 



XOR 
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where the "©" operation now indicates addition modulo 
3, can be decomposed into three operations, 
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The generalized Fourier transform, F, for one qutrit is 
defined by 



1 
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where j = 0, 1, 2. 

The phase gate for qubits, P qu b\j)m\k) n = 
exp(i7rjfc)|j) m |fc)„, is sandwiched by two Hadamard 
gates on the target qubit to give the controlled-NOT 
or XOR operation. The generalization of this gate to 
qutrits is P mn - It is completely specified by 



P mn \j)m\k) n = exp(2i7rj/c/3)|j) m |fc) r 
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We can generate this gate between 2 qutrits using a com- 
bination of single qutrit rotations and evolution under the 
two qutrit Hamiltonian H12 = 2-kJM\ ® M^- One pulse 
sequence for P mn is 
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where (xffi 



\l) tran- 



(^ n) ) = C^,(f,0) and (MM) e 



represents a 7r-pulse on the \k) 

sition of the nth ion i.e. 
represents a period of evolution under the Hamiltonian 
H12 for a time t such that 2ixJt — 9. While this is not 
a unique decomposition, the number of periods of evolu- 
tion under Hmm is minimal. The angles {a{\^ i=1 have 
been numerically determined and are shown in Tabled 



Refocussing techniques Q developed for NMR quan- 
tum computing are necessary here given the "always-on" 
nature of Hmm- For qubits the interaction is H$s — 
2ir J a z 2 \ The relations o a a z a a = — a z for a = x,y, 
where {<?i}i=x,y,z are the Pauli operators, are used to 
reverse the evolution under H$s 



e itH 3S = e - 



(17) 



Essentially the diagonal elements of a z are permuted by 
the two 7r-pulses. Combined with another period of evo- 
lution under H$s, the trace- less property of a z is ex- 
ploited so that nothing happens. In fact, all qubits cou- 
pled to the first qubit are refocussed by this. 

For qutrits the de-coupling procedure is basically the 
same. The spin-spin term Hmm = 2~kJ M\M 2 is com- 
posed of nine terms, 



MiM 2 = a H <8> M + (a 3 A 3 + a 8 A 8 ) ® M 



(18) 



By applying pulses which simultaneously permute the en- 
tries of A 3 , and A 8 , over three different permutations, in a 
manner similar to the qubit case, the combined evolution 
under the second term in l|18|) is removed. The sequence 
to refocus the MM evolution over an arbitrary dura- 
tion T, requires one to subdivide the duration into three 
equal-duration sub-evolutions (MM)g, (2wJT/3 = 6). 
One can then obtain, 
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and 
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The global phase factor, (j) — 3a^9, is irrelevant while 
the two single-qutrit pulses on the second qutrit are due 
to the first term in 1|18[) and can be easily reversed. The 
three periods of evolution required are expected since the 
matrices are 3-dimensional. 

Readout of the final state of the qutrit register takes 
two steps. The entire ion chain is illuminated with opti- 
cal radiation and the observed fluorescence spatially re- 
solved. The radiation frequency is chosen so that if the 
qutrits are projected onto |2) this is then detected. A 
7r-pulse is then applied on the |1) |2) transition of 
each ion and the ion string illuminated again. Fluores- 
cence now indicates projection onto |1) while its absence 
means the qutrit is in |0). 

We now an describe explicit example using the 171 Yb + 
ion. The hyperfine constant for m Yb + is A = 
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2ir 12.6GHz. In a magnetic field its levels are split as 
shown in Fig. |21 In fields of around 0.45T, the Rabi- 
Breit region, the good quantum numbers are F and Mp 
and our logical states are |0) = 1651 F = l,Mp = — 1), 

|1) = \6SiF = 1, Mp = 0), and |2) = \6SiF = 1,M F = 

2 2 

1). The transition frequencies loqi, and u>i2, are about 
3.7GHz and 8.9GHz, respectively while the differences 
between the resonance frequencies of two neighbouring 
ions in a trap of frequency 2tt 200kHz and a magnetic field 
gradient of b = lOOT/m are Sujqi s=s 11 MHz, and 6u>i2 ~ 
2MHz. The operator M n in © has an element which 
is approximately l/y/2. This approximation introduces 
a constraint on the size of the magnetic field gradient if 
we require M n to be constant over the ion chain within 
an accepted error cm- Numerical calculations in |27ll28|. 
give the minimum distance between ions in a chain as 

Az min (N) = 2.018 7 7V - 559 where 7 = (q 2 /4ne mu 2 )3 . 
As a conservative estimate, let us say that the ions are 

equally spaced at Sz = 1.5Az min (N). This implies that 

2 

b < O.OSTV -0 - 441 ^ 3 , limiting the size of the magnetic field 
gradient. On the other hand we need to be able to fre- 
quency discriminate between qutrit transitions on neigh- 
bouring ions, and be far enough away that no vibrational 
motion will be excited i.e. (dujQi/dz)5z < 2i//v + V\. Us- 
ing the numerical result 20] that = (2.7 + 0.57V)fi, 
this imposes a minimum size on the magnetic field gra- 
dient of b > 1.5 x lO" 9 ^ (3.27V 559 + 0.5A 1 - 559 ). For 



an axial trap frequency of v\ = 2ir 200kHz containing 10 
ions and cm = 0.01 the magnetic field gradient is limited 
by 30T/m < b < 200T/m. These limits also determine 
the maximum number of ions we can place in the trap, 
assuming other considerations related to the ratio of the 
axial and radial trapping frequencies are satisfied [23l |. 
For an axial trap frequency of v\ — 2n 200kHz, the max- 
imum number of ions is about 30 satisfying the above 
limits when the magnetic field gradient is 120T/m. The 
expression for the J-couplings in I© is the same as that 
in the qubit case. For 10 171 Yb + ions at a trap frequency 
of 2ir 200kHz and magnetic field gradient 120T/m gives 
a nearest-neighbour coupling of about J = 2?t 1.2kHz. 

In summary, we have presented a modification of previ- 
ous designs for ion quantum computation with magnetic 
field gradients but where the quantum information is now 
manipulated and stored in qutrits. A magnetic field gra- 
dient allows for individual qutrit addressing and intro- 
duces a qutrit-qutrit coupling for quantum logic. The 
scheme requires no additional physical resources. 
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TABLE I: The angles ai written in multiples of 7r required 
in order to execute the pulse sequence for the phase gate for 
qutrits given in 1171 . 
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FIG. 1: The qutrit levels {|0}, |1), |2}} under consideration 
with W12 and ojoi the resonant frequencies for the |1) <-> |2) 
and |0) «-» |1) transitions respectively. 




FIG. 2: The hyperfine Zeeman levels of an m Yb+ ion in a 
spatially varying magnetic field. 



